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ABSTRACT 

We analyse brightness profiles for 143 early-type galaxies in the Virgo and Fornax 
Clusters, observed with the Advanced Camera for Surveys on the Hubble Space Tele- 
scope. Sersic models are found to provide accurate representations of the global profiles 
with a notable exception: the observed profiles deviate systematically inside a character- 
istic "break" radius of Rb ~ 0.02^q'q^^ R^, where Re is the effective radius of the galaxy. 
The sense of the deviation is such that bright galaxies {Mb < —20) typically show 
central light deficits with respect to the inward extrapolation of the Sersic model, while 
the great majority of low- and intermediate-luminosity galaxies (—19.5 < Mb ^ —15) 
show central light excesses; galaxies of intermediate luminosities (—20 < Mb < —19.5) 
are generally well fitted by Sersic models over all radii. We show that the slope, 7', of 
the central surface brightness profiles, when measured at fixed fractions of Re, varies 
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smoothly as a function of galaxy luminosity in a manner that depends sensitively on 

the choice of measurement radius. We find no evidence for a core/power-law dichotomy, 
and show that a recent claim of strong bimodality in 7' is likely an artifact of the biased 
galaxy selection function used in that study. To provide a more robust characterization 
of the inner regions of galaxies, we introduce a parameter, A0.02 = log {Cg/ Cs) — where 
Cg and Cs are the integrated luminosities inside 0.02i?e of the observed profile and of 
the inward extrapolation of the outer Sersic model — to describe the central luminosity 
deficit (A0.02 < 0) or excess (A0.02 > 0). We find that A0.02 varies smoothly over the 
range of ~ 720 in luminosity spanned by the sample galaxies, with again no evidence 
for a dichotomy. We argue that the central light excesses in Mb ^ — 19 galaxies may 
be the analogs of the dense central cores predicted by some numerical simulations to 
form via gas inflows. 

Subject headings: galaxies: clusters: individual (Virgo, Fornax) -galaxies: elliptical and 
lenticular, cD-galaxies: nuclei: galaxies: structure 

1. Introduction 

Pioneering HST imaging studies of the centers of early-type galaxies suggested an apparently 
abrupt transition in central stellar density at Mb ~ —20.3 mag — the so-called "core/power-law 
dichotomy" (e.g., Ferrarese et al. 1994; Lauer et al. 1995). These findings prompted the widely held 
view that the bright ("core") and faint ("power-law") galaxies follow distinct evolutionary routes 
(e.g., Faber et al. 1997). However, the evidence for such a dichotomy has lessened (although not 
entirely disappeared) following more recent studies that identified a population of galaxies with 
intermediate properties (Rest et al. 2001; Ravindranath et al. 2001). The slope, 7', of the central 
surface brightness profile — usually parameterized as a "Nuker" law (essentially two power-laws that 
merge at a characteristic "break" radius; Lauer et al. 1995) and measured at the angular distance 
corresponding to the instrumental resolution — has traditionally been taken as a diagnostic of this 
behavior. 

However, using ACS profiles for 100 early-type galaxies belonging to the Virgo Cluster, Fer- 
rarese et al. (2006a) showed that Sersic models provide more accurate parameterizations of the 
global brightness profiles than do Nuker models (see also Graham 2004; Ferrarese et al. 2006c) 
and argued that the core/powcr-law dichotomy is an artifact introduced in part by the use of an 
inappropriate (i.e., power-law) parameterization of the outer profiles, combined with a tendency 
in previous work (which usually relied on HST brightness profiles of limited radial extent) to not 
properly account for the compact stellar nuclei found in low- and intermediate luminosity galaxies 
(e.g., Graham & Guzman 2003; Grant et al. 2005; Cote et al. 2006). 

In this paper, we use the best available imaging dataset — in terms of depth, radial coverage, 
angular resolution, completeness and homogeneity — to re-examine the central structure of early- 
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type galaxies. Our analysis relies on HST/ACS imaging for 100 early- type members of the Virgo 
(previously discussed in Ferrarese et al. 2006ab and Cote et al. 2006) and new HST/ACS imaging 
for 43 early- type members of the Fornax cluster (Jordan et al. 2007). Our principle finding is clear 
evidence for a continuous, systematic progression from central luminosity deficit {Mb < —20) to 
excess {Mb > —19) within a characteristic radius, approximately equal to 2% of the galaxy effective 
radius. We find no evidence for a "core/power- law" dichotomy. 

2. Observations 

HST images for 143 members of the Virgo and Fornax Clusters were acquired with the Ad- 
vanced Camera for Surveys (ACS, Ford et al. 1998) as part of the ACS Virgo (ACSVCS; GO-9401) 
and Fornax (ACSFCS; GO-10217) Cluster Surveys (Cote et al. 2004; Jordan et al. 2007). Surface 
brightness fluctuation distance measurements (Mci et al. 2005, 2007; Blakeslee et al. 2008, in prepa- 
ration) reveal the program galaxies to span luminosity ranges of ~ 545 (Virgo), ~ 345 (Fornax) 
and 720 (combined). All galaxies have early- type morphologies (i.e., E, SO, dE, dE,N or dSO) 
and are confirmed velocity members of their respective clusters. Images were taken in the Wide 
Field Channel (WFC) mode with a filter combination (F475W and F850LP) roughly equivalent to 
the g and z bands in the Sloan Digital Sky Survey (SDSS) photometric system. The images cover 
a roughly 200" x 200" field with 0"05 pixel^"*^ sampling and ~ O'.'l resolution. This resolution limit 
translates to physical scales of 8.0 and 9.5 pc for Virgo and Fornax, respectively. For the 21 Virgo 
galaxies brighter than Bt = 12 {Mb ~ —19.2), the profiles were extended in radius by matching the 
ACS profiles to those measured from SDSS g and z mosaics (Data Release 5; Adelman-McCarthy 
et al. 2007) generated using the procedures described in West et al. (2007). The ACSVCS sample 
is complete for early-type galaxies brighter than ~ 12 {Mb ~ —19.2) and 44% complete down 
to its limiting magnitude of S ~ 16 {Mb ~ —15.2). The ACSFCS sample, meanwhile, is complete 
down to its limiting magnitude oi B Ki 15.5 {Mb « —16.1) with the exception of a single galaxy 
(FCC161 = NGC1379) for which the data acquisition failed because of a shutter problem. 

Full details on the construction of the azimuthally averaged brightness profiles and the adopted 
fitting procedures (e.g., correction for dust obscuration, masking of background sources, the iden- 
tification of offset nuclei via centroid shifts, and the choice of weighting schemes and minimization 
routines, etc) are given in §§3.1, 3.2 and 3.3 of Ferrarese et al. (2006a) and §§3.1, 3.2, 4.1 and 4.2 
ofCoteetal. (2006). 

2.1. Parameterization of the Brightness Profiles 

Figures 1 and 2 shows g and z surface brightness profiles for nine representative galaxies 
belonging to each of the ACSVCS and ACSFCS samples. To parameterize the brightness profiles. 
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we begin by noting that Sersic (1968) models 

Is{R) = 4exp \-br^ 
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l/n 
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(1) 



generally provide accurate descriptions of the global profiles of the program galaxies — most im- 
portantly the downward curvature on large scales — with only three free parameters (/g, n and 
Re). However, the central regions of the profiles (typically i? < 100 pc for the brightest galaxies, 
and i? < 10 pc for the faintest), deviate significantly from the inward extrapolations of the Sersic 
profiles for both galaxies brighter than Mb — 20 (which show central light deficits) and for most 
galaxies fainter than Mb > —19.5 (which usually show central excesses). By contrast, most galaxies 
with —20 < Mb ^ —19.5 are reasonably well fitted by a single Sersic model over all radii, as previ- 
ously noted by Graham &; Guzman (2003). The discovery of a systematic progression from central 
light deficit to excess along the luminosity function was discussed extensively in Cote et al. (2006) 
and Ferrarese et al. (2006ab). This transition from deficit to excess at —20.0 < Mb < —19.5 
corresponds to Sersic indices of 3.9 < n < 3.5, based on the scaling relations presented in §4.3 of 
Ferrarese et al. (2006a; see their Eq. 27). 

To put this trend on a more quantitative footing, we parameterize the observed brightness 
profiles of the ACSVCS and ACSFCS galaxies in two different ways. First, given that the departures 
from Sersic models always occur on scales of a few arcseconds and usually much less — meaning 
that the innermost behavior is dominated by the instrumental point spread function (PSF) — 
a simple, PSF-convolved, power- law profile (with < 7' < 2) inside a break radius, usually 
provides an adequate representation of the observed profiles. The short dashed curves in each panel 
of Figures 1 and 2 shows such "core-Sersic" models (Graham et al. 2003), 
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fitted to the profiles in both bandpasses after convolution with the instrumental PSFs. Note that 
/' in Eq. (2) is related to the intensity. If,, at the break radius Rf, through the relation: 



I' = h2-^l'^ exp 



bn{2^/^Rb/Re 



l/n' 



(3) 



Thus, this choice of parameterization requires a total of five free parameters (/t, a, 7, n and Re). 
For comparison, the dotted curves in each panel of Figures 1 and 2 show the inward extrapolations 
of an outer Sersic component, illustrating the systematic evolution from deficit to excess. The 
arrows show the fitted values of Rf, for each galaxy. 

As pointed out above, the profiles in the innermost regions are dominated by the instrumental 
PSF, so a different parameterization of the inner component could provide equally acceptable fits, 
particuarly for the low- and intermediate-luminosity galaxies. Indeed, a detailed study of the 
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Fig. 1. — Surface brightness profiles for nine representative galaxies from the ACSVCS. Identifi- 
cation numbers from the Virgo Cluster Catalog (VCC) of Binggeli, Sandage & Tammann (1985) 
are given in the upper right corner of each panel. For each of these galaxies, which span a range 
of ~ 545 in blue luminosity, we show both the g and z profiles as the lower and upper points, 
respectively. The galaxies are ordered according to decreasing absolute blue magnitude which is 
recorded in each panel, along with their morphological types and the best-fit Sersic indices, n, for 
the galaxy measured in the g and z profiles (lower and upper labels). For each profile, we show the 
best-fit "core-Sersic" (short dashed curves) and "composite" model (long dashed curves; see §2.ip . 
The red and blue arrows show the break radii, R^, for the fitted core-Sersic models (in g and z, 
respectively) while the dotted vertical lines in each panel are drawn at 2% of the effective radius. 
The dotted curves show the inward extrapolations of the Sersic component that best fits the profile 
for R> Rf,. Note the smooth transition from luminosity "deficit" to "excess" as one moves down 

tVip InminnsitA;' fiinrtinn 
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Fig. 2. — Same as in the previous figure, except for nine representative galaxies from the ACSFCS. 
Identification numbers from the Fornax Cluster Catalog (FCC) of Ferguson (1989). The galaxies, 
which span a range of ~ 345 in blue luminosity, exhibit the same basic trends as those in the 
Virgo Cluster. Note the small break radius for FCC21 (= NGC1316 = Fornax A), a peculiar SO 
LINER/radio galaxy that shows clear evidence of recent merging (i.e., dust. Ha filaments, shells 
and ripples; Schweizer 1980). 
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nuclear brightness profiles of nearby galaxies (Ferrarese et al. 2008, in preparation) suggests that 
an excellent parameterization over all scales is given by a double Sersic model, 
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in which one component corresponds to the galaxy profile and the other to the central light excess 
(i.e., a compact spheroidal or flattened stellar component). In this case, a total of six parameters 
are needed to describe the entire brightness profile (/e,i, ni, Re,i, Ie,2, n2 and Re,2)- 

Thus, for our second ("composite") description of the ensemble profiles we adopt, on an 
object-by-object basis, one of three parameterizations that best fits the observed profile: (1) a 
core-Sersic model (Eq. 2); (2) a single, unbroken Sersic model (Eq. 1); or (3) a double Sersic model 
(Eq. 4). As mentioned above, the precise choice of parameterization is found to depend strongly on 
galaxy luminosity: all galaxies brighter than Mb ^ —20 are modeled with core-Sersic laws, while 
most galaxies slightly fainter than this (—20 < Mb ^ —19.5) are well fitted with a single Sersic 
model over all radii. For the majority of galaxies fainter than Mb ^ —19.5, a double Sersic model 
can accurately match the profiles on all scales. A small number of galaxies (< 10% of the sample) 
fainter than Mb ~ —17.5 are found not to require a second Sersic component (i.e., they contain 
no obvious central nucleus and their profiles are generally well fitted by a single Sersic model). We 
shall return to these interesting objects in ^ The long dashed curves in Figures 1 and 2 show 
these "composite models" , with the exact choice of parametrization listed in the lower left corner 
of each panel. The short and long dashed curves are generally indistinguishable, confirming that, 
for these particular galaxies, either approach yields an acceptable parameterization of the profiles 
on both large and small scales. Our conclusions are therefore robust to the choice of these two 
parameterizations. 



3. Results 

3.1. Central Surface Brightness Profile Slopes 

While it is well established that the brightest early-type galaxies differ from their fainter 
counterparts in terms of isophotal shape, ellipticity, kinematics and stellar populations (see, e.g., 
Bender et al. 1989; Kormendy & Djorgovski 1989; Caon et al. 1993; Ferrarese et al. 2006a; Emsellem 
et al. 2007) the strongest evidence for a bonafide dichotomy in their structural properties — rather 
than a continous variation along the luminosity function with some intrinsic scatter — has come 
from the behavior of their central brightness profile slopes (Ferrarese et al. 1994; Lauer et al. 1995; 
Rest et al. 2001; Ravindranath et al. 2001; Lauer et al. 2007). However, Ferrarese et al. (2006ac) 
noted that such slope measurements have often relied on heterogeneous archival HST images (i.e., 
different instruments, filters, and resolution limits) or on galaxy samples with ill-defined selection 
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functions and spanning wide ranges in distance (e.g., from 3.5 Mpc to ~ 320 Mpc in the case of 
the recent study by Lauer et al. 2007). If such a dichotomy has a physical origin, then it would 
seem more appropriate to measure the slope at either the same physical — rather than angular — 
radius, as was done in Ferrarese et al. (2006a), or at a point corresponding to a constant fraction 
of some characteristic scale radius in every galaxy (see also Ferrarese et al. 2006c). 

Panels (a-d) of Figure 3 illustrate this point by plotting the slope of the best-fit composite 
model at differing fractions of the effective radius of the galaxy (0.005, 0.01, 0.05 and 0.30) as a 
function of absolute blue magnitude0 That is to say, for each galaxy in our sample, we measure 
the instantaneous slope of the fitted models using the relations 

l/(on)-l 



(5) 

where the cS, S and dS subscripts refer to the core-Sersic, Sersic and double-Sersic parameterizations 
discussed in §2.1[ Note that this approach differs from the one adopted by Ferrarese et al. (2006a) 
who — recognizing that compact stellar mass concentrations become increasingly prominent among 
progressively fainter and lower surface brightness galaxies — preferred to measure the slope of the 
underlying galaxy component rather than using the combined (i.e., galaxy + nucleus) profiles. 
The photometric and structural scaling relations derived by Ferrarese et al. (2006a) corroborated 
previous reports of a continuum in galaxy properties fainter than Mb ~ —20 (e.g., see §4 and 5 of 
Ferrarese et al. 2006a and references therein), with no dichotomy between giants and dwarfs (see 
also the discussion in §1 of Graham & Guzman 2003). 

Panels (a-d) show the 'j'-Mb relations found when the profiles are parameterized as described 
in §2: a core-Sersic profile (Eq. 2) for the brightest galaxies, a double Sersic profile (Eq. 4) for most 
of the low- and intermediate-luminosity galaxies, and a Sersic profile in all other cases (mostly 
intermediate luminosity and a few low-luminosity galaxies). Four important points are illustrated 
in these panels. First, the overall trends defined by the Virgo and Fornax galaxies are indistin- 
guishable. Second, the behavior of the j'-Mb relatiorj^ is obviously sensitive to the exact choice 



^Absolute magnitudes for the ACSVCS galaxies are computed from the apparent magnitudes given in Cote 
et aL (2004), reddenings as described in Jordan et aL (2004) and SBF distances from Mei et al. (2007). For the 
11 Virgo galaxies lacking SBF distances, we assume (m — M) = 31.09 mag. For the ACSFCS galaxies, we use the 
apparent magnitudes and reddenings from Jordan et al. (2007) and SBF distances from Blakeslee et al. (2008, in 
preparation). Four ACSFCS galaxies without SBF distances are assigned the median Fornax distance. 

^We present these scaling relations in terms of Mb for two reasons. First, the scaling relations and dichotomies 
discussed in this paper have traditionally been expressed in this bandpass. Second, placing the ACSVCS and ACSFCS 
galaxies in the broader context of their cluster environments (i.e., see i]3.2|) is most straighforward using this bandpass 
since the wide-field, photographic (blue emulsion) surveys of Binggeli, Sandage & Tammann (1985) and Ferguson 
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Fig. 3. — (Panels a-d) Slope of the best-fit composite model, 7', measured at different fractions of 
Re (i.e., 0.005, 0.01, 0.05 and 0.30). Virgo and Fornax galaxies are shown by blue and red symbols, 
respectively. The behavior of the j'-Mb relation depends sensitively on the radius at which 7' is 
measured, but there is no evidence for a "core/power-law" dichotomy. The dashed curves show 
low-order Legendre polynomials that highlight the trends with magnitude. (Panel e) Dependence of 
^0.02 = ^og[Cg/ Cg) on galaxy magnitude for the composite fits (see text for details). Galaxies with 
central luminosity deficits have A0.02 < while those with central excesses have A0.02 > 0. There 
is a smooth transition from central deficit to excess with decreasing galaxy luminosity. A typical 
errorbar is shown on the left side of the panel. Open symbols denote those galaxies with dE/dIrr 
transition morphologies, dust, young stellar clusters and/or evidence of young stellar populations 
from blue integrated colors. 



- 10 - 



of measurement radius, particularly on the smallest angular scales. Third, the j'-Mb relation un- 
dergoes an unmistakable "inversion" between 0.01-0.05i?e: on the smallest scales, the profiles 
are found to steepen systematically as one moves down the luminosity function, whereas on larger 
scales, the profiles become progressively shallower with decreasing luminosity. And, finally, in no 
case is there clear evidence for a dichotomy or bimodality (see also Figures 1 and 2). 

3.2. Comparison to Lauer et al. (2007) 

The continuity of the ^'-M^ relations shown in the first four panels of Figure O while consistent 
with the finding of no core/power-law dichotomy by Ferrarese et al. (2006a), is apparently at odds 
with a recent detection of strong bimodality in 7' by Lauer et al. (2007). However, it is clear from 
this figure that the distribution of slopes found for any sample of galaxies will depend sensitively 
on both the choice of measurement radii and the luminosity distribution of the galaxies themselves. 

In the upper panel of Figure [U we plot as the open blue histogram the luminosity function, 
0gai, of the 143 galaxies that make up the ACS Virgo and Fornax Cluster Surveys. The single- and 
double-hatched histograms show the results for the ACSVCS and ACSFCS samples, respectively. 
The dashed red curve shows a Schechter function with a = —1.40 and B* = 9.8 {M^ = —21.4), 
the best-fit parameters for early-type (E+SO+dE+dSO) galaxies in the Virgo Cluster according 
to Sandage, Binggeli & Tammann (1985). The normalization of the Schechter function has been 
chosen to match the luminosity function of the ACSVCS and ACSFCS sample galaxies brighter 
than Mb ~ —19. Recall that both surveys are complete above this level and that the ACSFCS 
sample is complete for Mb ^ —16.1. 

The blue histogram in panel (h) shows the luminosity distribution of the 219 galaxies analysed 
by Lauer et al. (2007) i The Schechter function from panel (a)\s reproduced as the dashed red curve. 
The common element in the selection of these galaxies was the availability in the literature of Nuker 
model fits to brightness profiles derived from either WFPCl, WFPC2, NIC2 or NIC3 imaging. Two 
important properties of the Lauer et al. (2007) sample are worth noting. First, it is unrepresentative 
of the Schechter function form that provides a reasonable match to the luminosity distribution of 
the early-type galaxy populations in Virgo and Fornax in particular (Sandage et al. 1985; Ferguson 
& Sandage 1988) and, more generally, those of galaxies in both cluster or field enviroments (e.g., 
Schechter 1976; Loveday et al. 1992; Marzke et al. 1994; Blanton et al. 2003). Second, as Lauer 
et al. (2007) point out, their sample is bimodal in luminosity — a result that we confirm. A 



(1989) remain the most homogeneous photometric catalogs for Virgo and Fornax cluster galaxies. Homogeneous 
'curve-of-growth' photometry from wide-field ugrizJHK imaging for the ACSVCS galaxies will be presented in Chen 
et al. (2008, in preparation). 

^For comparison to the ACSVCS and ACSFCS samples, we have converted the Mv magnitudes given in Lauer 
et al. (2007) to Mb by assuming {B -V) ^ 0.96 for BCG/E galaxies and {B-V)^ 0.85 for SO galaxies (Fukugita, 
Shimasaku & Ichikawa 1995). A color of [B — V) = 0.85 was also adopted for the lone Sa galaxy in their sample. 
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KMM analysis (McLachlan & Basford 1988; Ashman, Bird & Zepf 1994) indicates a 0.3% chance 
that these data are drawn randomly from a unimodal Gaussian distribution (the null hypothesis). 
Assuming homoscedasticity and clipping seven galaxies with Mb > —17.5 (which serve to make 
the distribution even more bimodal), we find peaks at Mb = —21.68 and —19.30. 

What then is the expected distribution of 7' values for this dataset? As panels (a-d) of FigureE] 
demonstrate, the detailed form of the 7' distribution will depend not only on the program galaxy 
magnitudes, but also on the choice of radius, R/Re, at which 7' is measured. Lauer et al. (2007) 
choose to measure 7 at the data resolution limit (0'.'02, 0'.'04 or O'.'l depending on the instrument) 
for all galaxies in their sample, regardless of their distance or effective radius. According to the data 
in their Tables 1 and 2, R/Re thus varies between 0.00017 and 0.030 (a factor of 178) for the 150 
galaxies in their sample for which estimates of Re are available in their Table 1. To illustrate the 
behavior of the 7' distribution expected for their sample, we combine their Mb values, calculated as 
described above, with the "j'-Mb relation shown as the smooth polynomial in panel (a) oi Figure O 
i.e., for R/Re = 0.005, which falls within the range of R/Re values used by Lauer et al. (2007). To 
approximate intrinsic scatter and measurement uncertainties, we assume a fixed dispersion of ct(7') 
= 0.1 for Mb < —22, and rising linearly with galaxy magnitude to cr(7') = 0.5 at Mb = —16. 

The 7' distribution found in a typical simulation is shown by the heavy magenta histogram in 
panel (c). In this case, KMM identifies highly significant peaks at 7' ~ 0.17 and 0.70, with a less 
than 0.1% chance of drawing these data from a unimodal Gaussian distribution. Thus, the simulated 
distribution is found to be strongly bimodal despite the fact that the input ^' -Mb varies smoothly 
along the luminosity function. The peak values of the simulated distribution should be compared 
to the values of 0.11 and 0.70 that we find when applying KMM to the actual 7' measurements 
tabulated in Table 2 of Lauer et al. (2007), reproduced here as the thin blue histogram. Although 
the two distributions differ in the details (which is not surprising since, as mentioned above, Lauer 
et al. do not measure 7' at a fixed fraction of Re or any other fixed physical scale), our exercise 
suggests that the bimodality in 7' reported by Lauer et al. (2007) is likely an artifact of a sample 
defined by a bimodal luminosity distribution that is: (1) unrepresentative of real galaxy ensembles; 
and (2) strongly biased in favor of bright, but intrinsically rare, galaxies that are known to have 
nearly constant-density cores. For instance, 56 of the galaxies in the Lauer et al. compilation (26% 
of their sample) are brighter than VCC1226 (M49 = NGC4472), the brightest member of the Virgo 
Cluster. 



3.3. A Characteristic Radius: R}, 

To first order, Sersic models provide accurate representations of the global, azimuthally av- 
eraged brightness profiles for almost all galaxies in our sample, largely independent of luminosity, 
morphological type, prior classification as giant or dwarf, and the presence or absence of morpho- 
logical pecularities such as rings, shells and bars. However, the central deviations provide a notable 
exception to this rule. It is striking that, in all galaxies, the inner departure from the Sersic model 
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Fig. 4. (Panel a) Luminosity function, </)gai, for the combined sample of 143 early-type galaxies from the ACS Virgo 

and Fornax Cluster Surveys (histogram). The dashed red curve shows a Schechter function with parameters a = —1.40 and 
B* = 9.8 (Mg ~ —21.4), appropriate for E+SO+dE+dSO galaxies in the Virgo Cluster (Sandage, Binggeli & Tammann 1985). 
The normalization is chosen to match the program galaxy luminosity function at Mb ^ ^19 since both the Virgo and Fornax 
samples are complete above this luminosity. (Panel b) Luminosity function for the 219 galaxies in the compilation of Lauer 
et al. (2007). The sample is strongly non-representative of a Schechter function (the dashed curve from panel (a) is reproduced 
for comparison). A KMM analysis shows the sample to be bimodal in luminosity, with peaks at Mb = —21.68 and —19.30, and 
a 0.3% chance of drawing these data randomly from a unimodal Gaussian distribution (the null hypothesis). The dotted and 
solid curves show the separate Gaussians, and their sum, as determined by KMM. (Panel c) A comparison of the distribution 
of 7' values measured by Lauer et al. (2007) (thin blue histogram) with that expected (magnenta histogram) based on their 
program galaxy luminosities and the -j'-Mb relation measured at = 0.005-Re (see Figure 3, panel a). KMM reveals there to 
be a less than 0.1% chance of drawing these data randomly from a unimodal Gaussian distribution. 
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occurs at very nearly the same fraction of the effective radius. Using the core-Sersic parameteriza- 
tion (which can be apphed uniformly to all galaxies) we find the median logarithmic ratio of the 
break radius to the effective radius to be {log(Rh/Re)) = —1.68 with rms scatter 0.33. Thus, the 
inner deviations from the Sersic law are found to occur at a characteristic radius of Rb ~ 0.02Re, 
with some scatter (the it la range is 0.010-0.045) that does not appear to correlate with any galaxy 
property. As an illustration, the vertical dashed lines in each panel of Figures 1 and 2 are drawn 
at 0.02Re, where Re is measured independently in the g and z profiles. 

3.4. Light Deficits and Excesses: The A0.02 Parameter 

Because the departures from the inward extrapolation of the outer Sersic component occur at 
Rh ^ 0.02Re, we define a parameter, Ao.o25 as the logarithm of the ratio of the total luminosity 
of the best-fit galaxy model (core-Scrsic or composite), Cg, inside this radius to that of the outer 
Sersic component in this same region: A0.02 = ^og{Cg/Cs)- Note that this description of the 
central behavior of the brightness profiles has the advantage relative to 7' that it is an integral 
measurement, rather than a differential one, and is thus less susceptible to noise in the profiles, 
PSF uncertainties, and the extent to which the stellar nucleus may be resolved. Galaxies with 
central light deficits then have A0.02 < 0, those with excesses have A0.02 > 0, and those whose 
profiles are well represented by Sersic models over the full range in radius have A0.02 = 0. 

Panel (e) oi Figure 3 shows the dependence of A0.02 on absolute blue magnitude for the sample 
galaxies. The overall trends exhibited by the Virgo and Fornax galaxies arc the same: a gradual 
transition from light deficits (—0.5 < A0.02 ^ 0) for galaxies brighter than Mb ~ —20 to light 
excesses (0 < A0.02 ^ 1-2) for the faintest program galaxies. This description of the inner galaxy 
also allows a straightforward definition of nucleated vs. non-nucleated galaxies. For instance, 
adopting A0.02 > as the definition of a central stellar nucleus, one finds 69% of the combined 
sample galaxies, and 82% of all galaxies fainter than Mb ~ —19, to be nucleated. These estimates 
are consistent with the values for ACSVCS galaxies reported in Cote et al. (2006). 

We point out that the tendency for early-type galaxies of low- and intermediate luminosity 
to have central brightness profiles that fall above Sersic laws is not a new result: several previous 
investigators had noted that such behavior in isolated galaxies or using small datasets (e.g., Binggeli 
& Jerjen 1998; Kormendy 1999; Stiavelh et al. 2001; Graham & Guzman 2003). However, the 
systematic variation along the luminosity function from central light deficit to light excess — with 
no apparent structural dichotomies between "core" and "power-law" galaxies at Mb ~ —20, nor 
between giants and dwarfs at Mb ~ — 18 (c.f. Kormendy 1985) — has only been possible thanks 
to the uniformity, homogeneity and luminosity coverage of the ACSVCS and ACSFCS datasets. 
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Fig. 5. (Panel a) Comparison of WFPCl and ACS/WFC brightness profiles for VCC1440, an EO galaxy (Mb = -15.94) 

previously classified as having a "power-law" profile (Lauer et al. 1995; 2007). The upper datapoints show the deconvolved, major 
axis, F555W WFPCl profile of Lauer et al. (1995) and their fitted Nukcr model. The lower datapoints show the azimuthally 
averaged ACS F475W profile along with the PSF-convolved, double Sersic model (dashed and dotted curves) whose sum (solid 
red curve) best fits the observed profile. The datasets and models have been shifted apart by 0.75 mag for clarity. Residuals for 
both parameterizations are shown in the lower panel; note the S-shaped residuals for the Nuker model, indicative of a central 
nucleus (blue symbols; see also Ferrarese et al. 2006c). (Panel b) Predicted luminous mass profiles for the remnants of merging 
disk/halo galaxy models, showing the "dense stellar core" that forms as a result of gas dissipation and star formation. Panel 
(b) is adapted from the data/simulations shown in Figure la of Mihos & Hernquist (1994). (Panel c) Alternate presentation 
of the ACS brightness profile for VCC1440, showing the excess light inside 0.02_Re (A0.02 = 0.387 ± 0.006). Compare with 
panel (b). In panels (b) and (c), the abscissa spans the same range when normalized to (10"^ < R/Re < 10). 
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4. Discussion 

The central deficits in bright ellipticals have been traditionally interpreted as evidence for 
core evacuation by coalescing supermassive black hole (SBH) binaries (Ebisuzaki et al. 1991; Faber 
et al. 1997; Milosavljevic & Merritt 2001), although at least one recent study suggests that the 
dissipationless collapse of initially cold stellar distributions in existing dark matter haloes can also 
give rise to such features (Nipoti et al. 2006). The discovery of a systematic transition from 
central deficit to excess should provide an important new clue to the formation of both high- 
and low-luminosity early-type galaxies and the connection to SBHs. It is worth noting that the 
mass fraction contributed by the central nuclei in low- and intermediate-luminosity galaxies is 
indistinguishable (~ 0.2%) from that of SBHs in bright galaxies, suggestive of a shared formation 
path (Cote et al. 2006, Ferrarese et al. 2006ab; Wehner & Harris 2006; Rossa et al. 2006). On the 
other hand, McLaughlin & King (2006) point out that feedback by stellar winds and supernovae may 
also be capable of reproducing the observed scaling relations without a direct connection between 
the nuclei and SBHs. And, as Ferrarese et al. (2006b) note, there are at least some galaxies (such 
as M32) that contain both a central excess (Worthey 2004) and a SBH (Verolme et al. 2002). 

A lengthy discussion of the selection effects that prevented the ubiquitous nature of these 
nuclei from being fully appreciated with ground-based imaging was presented in Cote et al. (2006). 
Not surprisingly, the foremost advangtage of HST/ACS in this regard is its exceptional angular 
resolution. One might think, then, that the detection of nuclei in low- and intermediate-luminosity 
would have been straightforward in early high-resolution HST imaging, but this often proved not 
to be the case. For instance, the left panel of Figure [5] compares our ACS profile for VCC1440 — 
a low-luminosity elliptical previously classified as a "power-law" galaxy — to the WFPCl profile 
of Lauer et al. (1995, 2007)0 In this case, the limited radial extent of the WFPCl profile, and the 
adopted parameterization of the galaxy profile as a power-law, allowed the full extent of the central 
excess to go unrecognized, although its signature is evident in the S-shaped residuals of the fitted 
Nuker model (see also Ferrarese et al. 2006c). 

In hindsight, the realization that compact stellar mass concentrations reside at the centers of 
most low- and intermediate galaxies should perhaps have come as no surprise. Many theoretical 
studies have shown that such mass concentrations may be a generic outcome of the galaxy formation 
process. Proposed formation routes include repeated mergers of star clusters drawn to the galaxy 
center via dynamical friction, the growth of p oc r~'^l^ density cusps due to two-body relaxation 
of stars orbiting a central SBH (Bahcall & Wolf 1976) and/or central star formation driven by gas 
inflows — likely modulated by stellar evolution, mergers, IGM confinement, and feedback from 
stellar winds, supernovae or pre-existing AGN (e.g., Barnes & Hernquist 1991; Mihos & Hernquist 
1994; McLaughhn et al. 2006; Li et al. 2007; Hopkins et al. 2007). 



''Lauer et al. (2007) report 7' = 0.89 for this galaxy. We measure (7') = 1.08,1.40,0.79 and 1.31 at RjR^ = 
0.005, 0.01, 0.05 and 0.3, respectively. These radii are indicated by the tickmarks in Figure O 
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A comparison of the properties (e.g., sizes, luminosity functions, mass fractions, colors) of nuclei 
and star clusters in ACSVCS program galaxies argues against the first scenario as the dominant 
formation mechanism (e.g., see §5.2 of Cote et al. 2006). Likewise, Bahcall-Wolf "cusps" generated 
by SBHs do not seem viable since most of the nuclei are resolved and thus much more extended 
than the predicted cusps (e.g., Merritt &: Szell 2006). On the other hand, gas inflow models — 
either internally or externally modulated — seem more attractive as they have long predicted the 
formation of "dense stellar cores" in early-type galaxies. Consider the center and right panels 
of Figure [5] which compare the brightness profile of VCC1440 (the same galaxy shown in the left 
panel) to the simulations of Mihos & Hernquist (1994). These authors showed that compact central 
starbursts are found in the remnants of disk/halo mergers and went on to argue that this central 
component would lead to "a break in the mass profile at small radii (~ 2%i?e)"- As Figure[5]and the 
discussion in ^3.31 show, this may turn out to be a remarkably prescient prediction. The fact that 
the deviations from the fitted Sersic models occur on roughly the same (~ 2%i?e) scale for both 
the central excesses and central deficits (which are usually attributed to binary SBH evolution; 
see above) presents something of a puzzle since it is unclear why these two unrelated physical 
mechanisms would give rise to deviations on the same (fractional) scale. 

In any case, it will be important to test the gas inflow scenario envisaged by Mihos & Hern- 
quist (1994) by measuring star formation and chemical enrichment histories for the central regions 
of bright (deficit) and faint (excess) galaxies alike. Both the measured color profiles (Ferrarese 
et al. 2006a) and the model-dependent integrated colors (Cote et al. 2006) of the nuclei and galax- 
ies in the ACSVCS suggests that stellar population differences between the two components may 
indeed exist, with the nuclei in the faintest systems generally being bluer than the underlying galaxy. 
On the other hand, at least some intermediate-luminosity galaxies show evidence for surprisingly 
red central components (e.g., g — z > 1.5). 

Finally, we point out that the open symbols in panel (e) oi Figure [3] denote galaxies with 
"dE/dIrr transition" morphologies: i.e., dust, young stellar clusters and/or evidence of young 
stellar populations from blue integrated colors (Ferrarese et al. 2006a; Cote et al. 2006). Such 
galaxies typically lack conspicuous central nuclei at their photocenters. The absence of a central 
excess may suggest that the process of nucleus-building — perhaps through delayed gas inflows and 
subsequent star formation — has yet to occur in these apparently young galaxies. 
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